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Abstract 

Recent results from CLEO on B — > Ktt indicate that the phase 7 may be 

substantially different from that obtained from other fit to the KM matrix 

elements in the Standard Model. We show that 7 extracted using B — ► Ktt, tttt 

is sensitive to new physics occurring at loop level. It provides a powerful 

method to probe new physics in electroweak penguin interactions. Using 

effects due to anomalous gauge couplings as an example, we show that within 

the allowed ranges for these couplings information about 7 obtained from 

B — ► Ktt, tttt can be very different from the Standard Model prediction. 
PACS numbers:11.30. Er, 12.60. Cn, 13.20. He, 13.25. Hw 
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The CLEO collaboration has recently measured the four B — > Kir branching ratios with 
0, Br{B ± -> n^K ) = (1.82t°;f ±0.16) x lO" 5 , Br(B ± -> i^vr ) = (1.21±g;|g±8^) x 10~ 5 , 
Br(B -> K±7tT) = (1.88±g:|i ± 0.13) x 10~ 5 and 5r(5 -> #°7r°) = (1.48±g:|?±8;it) x lO" 5 . 
It is suppressing that these branching ratios turn out to be close to each other because naive 
expectation of strong penguin dominance would give R = Br(K ± ir°) / Br{K'K ± ) ~ 1/2 and 
model calculations for Br(B° — > i^°7r°) would obtain a much smaller value. The closeness 
of the branching ratios with charged mesons in the final states may be an indication of large 
interference effects of tree, strong and electroweak penguin interactions f2|. It has been shown 
that using information from these decays and B^ — > 7r ± 7r° decays, the phase angle 7 of the 
KM matrix can be constrained || and determined in the Standard Model (SM). Using 
the present central values for these branching ratios we find that the constraint obtained on 
7 may have potential problem with 7 = (59. 5ltf ;|)° obtained from other constraints ||. 

If there is new physics beyond the SM the situation becomes complicated It is not 
possible to isolate different new physics sources in the most general case. However, one can 
extract important information for the class of models where significant new physics effects 
only show up at loop levels for B decays |||7|]. In this paper we study how new physics of 
the type described above can affect the results using anomalous three gauge boson couplings 
as an example for illustration. 

New physics due to anomalous three gauge boson couplings is a perfect example of models 
where new physics effects only appear at loop level in B decays. Effects due to anomalous 
couplings do not appear at tree level for B decays to the lowest order, and they do not affect 
CP violation and mixings in K° — K° and B — B systems at one loop level. Therefore they 
do not affect the fitting to the KM parameters obtained in Ref. 0. However they affect 
the constraint on and determination of 7 using experimental results from B — ► Ktt, titt, and 
affect the B decay branching ratios. 

The effective Hamiltonian H eff = {G F /^/2)[V* h V uq {c x Ox + c 2 2 ) - V£V tq Ei= 3 -io cA] 
responsible for B decays have been studied by many authors We will use the values of 
Ci obtained for the SM in Ref. |§ with 
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ci = -0.313, c 2 = 1.150, c 3 = 0.017, c 4 = -0.037, c 5 = 0.010, c 6 = -0.046, 

c 7 = -0.001a em , c 8 = 0.049a em , c 9 = -1.321a em , c w = 0.267a em . (1) 

The Wilson Coefficients are modified when anomalous couplings are included. They will 
generate non-zero 03-10 |1U|| . Their effects on B — > Ktt mainly come from C7_io- To the 
leading order in QCD corrections, the new contributions to cf_ 10 due to various anomalous 
couplings are given by, 



C7 C /a em = -0.287Ak 7 - 0.045A 7 + 1.397A#f - 0.145#f , 
C8 C /a em = -0.082Ak 7 - 0.013A 7 + 0.391A#f - 0.041#f , 
Cg c /a em = -0.337Ak 7 - 0.053A 7 - 5.65lA#f + 0.586#f , 
c^ c /a em = 0.069A/t 7 + 0.011A 7 + 1.143A#f - 0.119#f . (2) 

In the above we have used a cut off A = 1 TeV for terms proportional to Ak 7 and Agf . 
Contributions from other anomalous couplings are suppressed by additional factors of order 
(ml, m'^)/m^ v which can be safely neglected. There are constraints on the anomalous gauge 
couplings LEP experiments obtain [0, -0.217 < Ak 7 < 0.223, -0.158 < A 7 < 

0.074, and -0.113 < Agf < 0.126 at the 95% confidence level. Assuming that #f is the 
same order as Agf , it is clear that the largest possible contribution may come from Agf . 
In our later discussions we consider Agf effects only. 

To see possible deviations from the SM predictions for B — > data, we carried out a 
calculation using factorization following Ref. |0 with V us = 0.2196, V c b ~ — V ts = 0.0395, 



\V u b\ = 0.08V c b |T3 |. The branching ratios as functions of 7 are shown in Fig. 1. In this 
figure we used m s = 100 MeV which is at the middle of the range from lattice calculations 
T5| and the number of colors A^ c = 3. Since penguin dominates the branching ratio for 



Br(B + — > K°7T + ) which is insensitive to 7, we normalize the branching ratios to Br(B + — > 
K°ii + ) to reduce possible uncertainties in the overall normalization of form factors involved. 

From Fig. 1, we see that the central values for the branching ratios for the ones with at 
least one charged meson in the final states require the angle 7 to be within 75° to 80° rather 



than the best fit value 'jbest = 59.5° in Ref. ||. Larger 7 is also indicated by other rare B 
decay data ||16|| . When effects due to Agf is included the situation can be relaxed. The 
effects of Agf on B + — > K°n + , K + tt~ are very small, but are significant for B + — > K + n° 
and B — > K°7r°. With positive Agf in its allowed range, it is possible for the relative 
ratios of Br(B + — > K + tt°) to the other charged modes to be in agreement with data for 
7 = lbest . We note that Agf does not affect the ratio Br(B° -> K+n-)/Br(B+ -> K°n+). 
Its experimental value prefers 7 to be close to 75°. Of course we also note that this situation 
can be improved by treating N c as a free parameter to take into account certain non- 
factorizable effects. We find that with N c w 1.35, the central experimental values for the 
branching ratios of B decays into charged mesons in the final states can be reproduced for 
7 = 7&esi- It is not possible to bring Br(B° — > K°tt°) up to the experimental central value 
even with allowed Agf and reasonable values for N c and m s . 

If the present experimental central values will persist and factorization approximation 
with N c = 3 is valid, new physics may be needed. Needless to say that we have to wait 
for more accurate data to draw firmer conclusions. Also due to our inability to reliably 
calculate the hadronic matrix elements, one should be careful in drawing conclusions with 
factorization calculations. However, we would like to point out that data on rare B to Kir 
decays may provide a window to look for new physics beyond the SM. Of course, to have a 
better understanding of the situation one needs to find methods which are able to extract 
7 in a model independent way and in the presence of new physics. In the following we will 
analyze some rare B to Kn decay data in a more model independent way. 

Model independent constraint on 7 can be obtained using branching ratios for B ± — > Kn 
and B^ — > -kit from symmetry considerations. This method would only need information 
from charged B decays to Kn and tttt, and therefore this method is not affected by the 
uncertainties associated with neutral B decays to Kn modes. Using SU(3) relation and 
factorization estimate for SU(3) breaking effect among B ± -> Kn and B ± — > 7T7T, one 
obtains 
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A(n+K°) + V2A(n°K + ) = eA(n + K°)e iA *^-^, 

r 3 [Vefe|[Vea| Cg + C W , _ 3 [Vf&|[Vjd| ia Cg + C 10 

2 |Kfe||K s | ci + c 2 2|Kb||v4rf| ci + c 2 

_ ^\Vus\f K \A(^rr°)\ 

~ \V ud \ fn \A{ir+K°)y {6) 
where A0 is the difference of the final state elastic re-scattering phases 4*3/2,1/2 for I — 
3/2, 1/2 amplitudes. For f K /f n = 1.22 and Br(B ± -> vr ± 7r°) = (0.54±g;|J ±0.15) x 10~ 5 , we 
obtain e = 0.21 ± 0.06. 

The parameter S' EW is of order cg^o/ci^ which is much smaller than one and will be 
neglected in our later discussions. Sew is a t rue measure of electroweak penguin interactions 
in hadronic B decays and provides an easier probe of these interactions compared with other 
methods [[I?]. In the above contributions from 07$ have been neglected which is safe in the 
SM because they are small. With anomalous couplings this is still a good approximation. In 
general the contributions from cj^ may be substantial. In that case the expression becomes 
more complicated. But one can always absorb the contribution into an effective Sew- I n the 



SM for r v = \V ub \/\V cb \ = 0.08 and = 0.2196 0, S EW = 0.81. Smaller r v implies larger 
Sew- Had we used r v = 0.1, Sew would be 0.65 as in Ref. With anomalous couplings, 

we find 

6 EW = 0.81(1 + 0.26Ak 7 + 0.04A 7 + 4.33A#f - 0.45#f). (4) 

The value for Sew can be different from the SM prediction. It is most sensitive to Agf . 
Within the allowed range of —0.113 < Agf < 0.126 ]l^], Sew can vary in the range 0.40 ~ 
1.25. 

Neglecting small tree contribution to B + — > tt + K°, one obtains 

K + r 2 )/2-l-6 2 (l-^) 

cos7 = S EW tti a ± , 1 — S ' 5 

2e(cos A(p + eS E w) 

— r 2 = 4esin A0sin7, (6) 
where r 2 = ABr^K^ /[Br(7T + K°) + Br(n'K )] = 1.33 ± 0.45. 



If SU(3) breaking effect is indeed represented by the last equation in (|3|), and tree con- 
tribution to — > n^K is small, information about 7 and Sew obtained are free from 
uncertainties associated with hadronic matrix elements. Possible SU(3) breaking effects 
have been estimated and shown to be small P,f4],|T8| . The smallness of the tree contribution 



to B ± — > Ktt^ is true in factorization approximation and can be checked experimentally 
19fl . The above equations can be tested in the future. We will assume the validity of Eq. 



(|) and study how information on 7 obtained from B — > Kir, im decays depends on Sew- 

The relation between 7 and Sew is complicated. However it is interesting to note that 
even in the most general case, bound on cos 7 can be obtained. For A = (r\ + r 2 _)/2 — I — 
e 2 (l — S 2 EW ) > 0, we have 

cos 7 < 5 EW - 9 cf1 \ — c, or C0S7 > S EW - , - -= — r- (7) 
2e(l + eS E w) 2e(-l + eS E w) 

The sign of A depends on r±, e and Sew- As long as r± > 1.07, A is larger than zero at the 

90% CL. in the allowed range for e and any value for Sew- F° r smaller r±, A can change sign 

depending on Sew- F° r A < 0, the bounds are given by replacing <, > by >, < in the above 

equations, respectively. We remark that if r\ < 1, one can also use the method in Ref. [2C] 

to constrain 7. The above bounds become exact solutions for cos A0 = 1 and cos A0 = — 1, 

respectively. For e << 1, one obtains the bound | 00S7 — Sew\ > ( r + + r 2 — 2)/(4e) in Ref. 

I- 

We will use the central value for e and vary r\ in our numerical analysis to illustrate how 
information on 7 and its dependence on new physics through Sew can be obtained. The 
bounds on cos 7 are shown in Fig. 2 by the solid curves for three representative cases: a) 
Central values for e and r\\ b) Central values for e and la upper bound r\ = 1.78; and c) 
Central value for e and la lower bound r\ = 0.88. For cases a) and b) A > 0, and for case 
c) A < 0. 

The bounds with | COS7I < 1 for a), b) and c) are indicated by the curves (al, a2), (b) 
and (cl, c2), respectively. For ) and c) there are two allowed regions, the regions 

below (al, cl) and the regions above (a2, c2). For case b) the allowed range is below (b). 



When r± decreases from la upper bound to la lower bound, one of the boundaries goes up 
from (b) to (al) then moves to (c2). And the other boundary for case b) which is outside 
the range moves to (a2) then goes down to (cl). In case a), for Sew = 0.81(0.65) we find 
cos 7 < 0.18(0.015) which is way below the value corresponding to cos 7^4 ~ 0.5. With 
Agf = 0.126, cos 7 can be close to 0.5. For larger r± the situation is worse. This can be 
seen from the curves for case b) where cos 7 < for Sew up to 1.5. For smaller r\ the 
situation is better as can be seen from case c). In this case there are larger allowed ranges. 
7 ~ Ibest can be accommodated even by the SM. 

When the decay amplitudes for — > Kit, — > 7r ± 7r and the rate asymmetries for 
these decays are determined to a good accuracy, 7 can be determined using Eq. ([|) and its 
conjugated form. The original method |^TJ using similar equations without the correction 



Sew from electroweak penguin is problematic because the correction is large [22]. Many 



variations involving other decay modes have been proposed to remove electroweak penguin 



effects P3| . Recently it was realized || that the difficulties associated with electroweak 
penguin interaction can be calculated in terms of the quantity Sew- 

This method again crucially depends on the value of Sew- The solution of cos 7 corre- 
sponds to the solution of a fourth order polynomial in cos 7. Using Eqs. (|5|) and <^), we 
have 

(1 - cos 2 7) [1 - (ttj^ r " «W) 2 ] - ^ '/P 2 = 0. (8) 

2e(S E w - cos 7) 16e^ 

The solutions depend on the values of r± and e which are not determined with sufficient 
accuracy at present. To have some idea about the details, we analyze the solutions of cos 7 as 
a function of Sew for the three cases discussed earlier with a given value for the asymmetry 
A aS y — ( r + — r - ) / ( r + + r2 _ ) • in Fig. 2 we show the solutions with A asy = 15% for illustration. 
The actual value to be used for practical analysis has to be determined by experiments. The 
solutions for the three cases a), b) and c) are indicated by the dashed, dot-dashed and dotted 
curves in Fig. 2. In general there are four solutions, but not all of them are physical ones 
satisfying | 00S7I < 1. 



For case a), two solutions are allowed with Sew- To have cos 7 > Sew has to be larger 
than 0.7. Whereas cos 7 ~ cos7b est would require Sew to be larger than 1.2 which can not 
be reached in the SM, but is possible for non-zero Agf in its allowed range. With smaller 
r±, cos 7 > can be solution with smaller 5ew and can even have cos 7 = cos 7;^. This can 
be seen from the dotted curves in Fig. 2 for case c). For larger r\ in order to have solutions, 
larger 5ew is required. For case b) Sew must be larger than 1.4 in order to have solutions. 
These regions can not be reached by SM, nor by the model with Agf in the allowed range. 

We also analyzed how the solutions change with the asymmetry A asy . With small A asy , 
the solutions are close to the bounds. When A asy increases, the solutions move away from 
the bounds. The solutions below the bounds (al), (b) and (cl) shift towards the right, and 
the bounds (a2) and (c2) move towards the left. In all cases the solutions with | COS7I 1 
become more sensitive to Sew an d | COS7I becomes smaller as A asy increases. In each case 
discussed the solutions, except the ones close to | C0S7I = 1, in models with non-zero Agf 
can be very different from those in the SM. It is clear that important information about 7 
and about new physics contribution to Sew can be obtained from B ± —> Kir, irir decays. 

We conclude that the branching ratios of B — > Kir decays are sensitive to new physics 
at loop level. The bound on 7, extracted using the central branching ratios for B ± — > Kir 
and information from B^ — > ir^ir, is different from that obtained from other experimental 
data. New physics, such as anomalous gauge couplings, can improve the situation. Similar 
analysis can be applied to any other model where new physics contribute to electroweak 
penguin interactions. The decay modes, B^ — > Kir, irir will be measured at various B 
factories with improved error bars. The Standard Model and models beyond will be tested 
in the near future. 
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FIG. 1. CP-averaged branching ratios normalized to Br(B + — » K°tt + ) vs. 7 for B + — > K°Tt + 
(dashed), B+ -> K+ir (solid), 5° -> K + tt- (dot-dashed), and B° -> i^°7r° (dotted) for the 
SM with m s = 100 MeV. The curves 01, b\ and 02, &2 are for Agf equal to —0.113 and 0.126, 
respectively. 




FIG. 2. Bounds on and (solutions for) cos 7 vs. 5ew- The curves al, a2 (dashed), b 
(dot-dashed), and cl,c2 (dotted) are the bounds (solutions) on (for) cos 7 as functions of Sew 
for the three cases a), b) and c) described in the text, respectively. 



11 



